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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1286 


METHOD OF SUCCESSIVE APPROXIMATIONS FOR THE SOLUTION OF 
CERTAIN PROBLEMS IN AERODYNAMICS* 

By M. E. Shvets 


The approximate solution of certain problems in the boundary-layer 
theory is considered herein. The method used is a combination of 
boundary-layer and successive-approximation methods . The essential 
character of the method will be evident from the problems considered. 

1. Solution of diffusion equation. - A number of physical processes 
are described by the diffusion equation 


u . f — a- 

Sx’ 5 Z »2 


and the boundary conditions 

q' (0,x) = Q = constant 

q' (°°,x) = 0 

q' (z,0) = 0 


( 1 . 1 ) 


( 1 . 2 ) 


where q is the concentration and u the velocity. Expressed in non- 
dimensional magnitudes 

q' = Qq 

i k' 
z ' = — z 
u' 


x» = 


U’ 


*"0 Pribl i zhennnm Reshenl i Nekotnrykh Zadach Ci d-rndi nam-i Png-rnrr! nViri<igO 
Sloya." Prikladnaya Matematika i Mekhanika, Vol. 13, No. 3, 1949, 
pp. 257-266. 
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equation (1.1) assumes the form 



(1.3) 


The solution for boundary conditions (equation (1.2)) is 


q = 1 




(1.4) 


where ^ is the error function. 


This solution gives the following expression for the diffusion 
velocity : 



(1.5) 


In considering the approximate solution of this problem, function 
6 termed the "thickness of the boundary layer" is introduced and the 
condition for z = o° is written as the condition for z = 6. The fol- 
lowing conditions are thus obtained? 

q (0,x) = l\ 

> ( 1.1 

q (S,x) = 0 J 

The thickness of the boundary layer 5 is an unknown function of 
the variable x; in order to determine it, the condition 


5 

must be assumed. 

With the aid of successive approximations, an attempt will be made 
to obtain q: 
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'N 


9 (1) - 4 


0 


<1 (2) = 4 0 + ll 
<l v ' = 4 q + 1l + 12 • • • 


q.( n+1 ) = q. 0 + 


■ • + 9n J 

where each of the functions is determined by equations (1.9) 


a 2 


<10 


az i 


= 0 




a 2 


*1 


a z ‘ 


o 2 dc 


a z ‘ 


= a, 


5 a n + 

a z 2 


0 
k 
% 


J 


( 1 . 8 ) 


(1.9) 


where the dot denotes differentiation with respect to x and satisfies 
the boundary conditions of equations (1.6). 

The operation of finding q from, the given £ with the aid of 
equations (1.9) and all boundary conditions of the problem are denoted 

®(t). 

Equations (1.9) then assume the form 

■\ 

, z 

*0 " 1 “ & 

a x = $(4 0 ) . • • 

Vi = 


(1.10) 
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The process of successive approximation considered gives the series 

q = q. 0 + <i>(4 0 ) + $(4!) + . . . + ^(qj (1.11) 


Stopping at the second approximation gives the following expression 
for q: 


, z zS 
1 '5 + -6 


Hi - ') 


( 1 . 12 ) 


For boundary- layer thickness from the condition of equation (1.7) 

55 = 3 
5 = 2.45 


(1.13) 


For the value of 6 determined, the following equations are 
obtained from equation (1.12) for the diffusion velocity and the con- 
centration q: 

'\ 

_ 1 + 5 ~ 0.61 
0~ 6 6 ~ 



•, £ 2.45 , 

4 " 1 " 2.45 + 12 ^ 


6 " 1 


( 2 * * 5 ‘ yr) 


(1.14) 




By comparing the diffusion velocity computed by the exact equation 
(1.4) and the approximate equation (1.14) with the corresponding curves 
in figure 1, the effectiveness of the approximation method can be seen. 

2. Approximate integration of boundary -layer equation of flat 

plate in incompressible flow. - A flat plate moves with constant velocity 

U. The equations of motion and continuity for the laminar boundary 

layer will be 
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5. 


, , ciu* . Su' 

u 5F- + t SF = u 


By ,£: 




du' , hr' _ n 


(2.1) 


J 


where V is the , kinematic viscosity. 

Elimination of the vertical component of velocity gives 


u 


5 2 u* 

Sy.2 


= U' 


5u* 

§x' 


5u 

3F1 


f 


Bu' 

SF dy 


( 2 . 2 ) 


The nondimens ional magnitudes 

x’ = Zx 

V' = V 


'S 


U 

^/E 


y ~ YK y 


u' = Uu 


(2.3) 


J 


are introducted where Z is the length of the plate and R is the 
Reynolds number. Equation (2.2) assumes the form 


d 2 u _ Bu 
^2 ‘ 


Bu 

5^j 


'y 


Bu 

5^ 


dy 


(2.4) 


The solution of equation (2.4) must satisfy the conditions 

u = 0 for y = 0 
u = 1 for y = 6 


(2.5) 


The first approximation is obtained by solving the equation 

P g 

B u/By = 0 for the boundary conditions of equations (2.5). This 
approximation has the form u = y/s. By substituting the first 
approximation for u in the right side of equation (2.4) and carrying 
out two quadratures, the second approximation is obtained! 


k. ( i _ z_\ 

5 -24 I 1 6 3 J 


vl = L + %L 

— r* ” a 


( 2 . 6 ) 
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As is shown by computation, the third approximation is not required. 
The magnitude 5 is determined from the condition (du/dy)g = 0. 

S = 4 -fx" (2.7) 


The frictional stress T according to equation (2.6) is 

_/c3u\ _ _1 _S 0.555 

T 'H ~ s + 24 ~ yr- 


( 2 . 8 ) 


or in dimensional magnitudes 

T’ 



(2.9) 


This equation differs from the exact numerical coefficient, which 
according to Blasius (reference 1) is 0.552. 


Substituting the value of 5 in equation (2.6) for the velocity 
profile yields 



( 2 . 10 ) 


The vertical component of the velocity is determined from the con- 
dition of continuity 



( 2 . 11 ) 


In figure 2 are shown the velocity profiles obtained by the approx- 
imate and exact solutions. The exact solution is mathematically rather 
complicated. 


5. Motion of fluid in laminar boundary layer with longitudinal 
pressure drop. - In this case, the motion is described by the equation 


d 2 u Su 

3^2 = U ST 



■\r 

0 


Bu . 


dx 


with the boundary conditions 

u = 0 for y = 0 


u = U for y = 5 
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If the calculation is restricted to the second approximation, 

= g? s 2 u 4 - 12 1 2 + no - ^ u 4 - e) + 1 u = y/eo 


The equation for determining the houndary -layer thickness, which is 
obtained from the condition (Su/S ^ = 0 has the form 

| 5 2 fr + - US5 = 1 

If it is assumed that 6 =0 for x = 0 

'x 


- — r 

D 6 Jo 


l^dx 


For the frictional stress at the wall 


- f I 1 + 4 D5ai )’ 1/2 f + C H 


The point of separation of the boundary layer is obtained from the 
condition 

4U f U^dx = - 1 


U 6 v^ 


0 


4. Heat— transfer in laminar boundary layer on plate in incompres- 
sible fluid, constants of which are independent of temperature. - If 
the addition of heat due to friction is neglected and the process is 
considered steady, the equation of heat conductivity is obtained: 


u* 


St 

Sx' 


+ v* 


St _ v a 2 T 

“ By ' 2 


(4.1) 


When the nondimens ional magnitudes (equations (2.3)) and the non- 
dimensional temperature 



(4.2) 
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where is the absolute temperature of the plate, and the abso- 

lute temperature of the fluid on the edge of the boundary layer, are 
introduced equation (4.1) assumes the form 


u 


5 


+ V 




Sy 


(4.3) 


where P is the Prandtl number. The boundary conditions will be 


II 

o 

for 

o 

ii 

a = i 

for 

y = 5 


If the values of the velocities according to Blasius (reference l) 
are substituted in equation (4.3), the equation is solved and leads to 
the following result : 


Hi) 



(4.4) 


where 


C(P) 



(4.5) 


As is seen from these relations, for a given Prandtl number P the 
temperature and the flow of heat can be obtained only by the method of 
numerical integration. The solution is rendered considerably more com- 
plicated if, in equation (4.1), the effect of the heat due to the work 
of friction is considered. 


By the interpolation of Pohlhausen in the interval 0.6 < P < 15.0, 
the function C(P) may be represented In the form 

C(P)= 0.332^P" 
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The problem of the heat -transfer of the plate has been the subject 
of many investigations, for example references 2 and 3. Piercy and 
Schmidt (reference 4) applying an approximate method similar to the 
method of Oseen arrived at the result 

C(P) = 0.404 W 

This result is about 20 percent greater than the result of 
Pohlhausen. 

The problem under consideration will now be solved by the approx- 
imate method. As a first approximation, d = y/g. 

Substituting the values of the component velocities obtained in 
section 2 and the derivatives of the temperature computed by the first 
approximation in the left side of equation (4.3) yields the second 
approximation : 


= 





(4.6) 


where 



From the condition (c)T/Sy)g = 0, boundary -layer thickness 
obtained 



(4.7) 


is 


(4.8) 


At a Prandtl number P = 1, the boundary-layer thickness is 
10 percent less than in the dynamic problem of the flow about a plate. 
This difference is explained by the fact that the solution (4.6) is 
approximate. If P = 0.8 (air) is assumed, the boundary-layer thickness 
will be 


5 = 4 -yfT 

By differentiating equation (4.6) with respect to y and setting 
y = 0, the following expression for the flow is obtained: 



II II III 


I 1 1 


■ ■■■■III I^HI 
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( St\ 3 l 

4 -Wo ■ °- S68 

Wien this result is compared with the results obtained by others 
(Piercy and Schmidt, reference 4), it is found to be closer to the 
accurate value. 

Substituting equation (4.8) in equation (4.6) yields the tempera- 
ture profile 



In figure 3 are given the profiles of the nondimens ional tempera- 
ture *9 for the values P = 1 and P = 7 . The solid curves are com- 
puted by the approximate formula (4.6) and the dashed ones by the 
accurate formula. These curves lie close to each other. 

5. Cooling of heated sphere by laminar flow of fluid at small 
values of Reynolds number. - The problem of the slow steady laminar flow 
about a sphere was considered by Stokes (reference 5). The magnitude of 
the tangential component of the velocity of the fluid at the distance r 
from the center of a spherical drop of radius r Q is determined by the 
expression 

u 0 = U sin 0 £ + | | ^ (5.1) 

The magnitude y is introduced by the relation r = r n + y. 

For small value of y (near the surface of the sphere) 0 

u 0 = 3U f sin 9 (5 . 2 ) 


In accordance with Leibenson (reference 6), the equation of heat 
conductivity is presented in the form 


„ St _ v 3 2 t 

S7 


(5.4) 


or, substituting the value of Uq 

U 


5 sin 9 v BT 5 T 

if 


kd 


(5.5) 
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The solution of this equation must satisfy the following boundary 
and initial conditions ; 



T = 0 


T = 0 


for 

for 

for 



(5.6) 


For simplification of the integration of equation (5.5), in accor- 
dance with Leibenson (reference 6), the magnitude (3/kd) sin 6 is at 
first assumed constant. The following variables are introduced: 

y = 2 

kd x 

U 6 sin 0 

Instead of equation (5.5) equation (5.8) is obtained; 



St 

51 



St 

3f 


(5.8) 


The solution of this equation is now obtained by the approximate 
method. As previously, the boundary-layer thickness is denoted by 5. 
If the calculation is restricted to the second approximation, the 
approximate solution of equation (5.8) is obtained; 


T = 1 - 





(5.9) 


Determining the boundary -layer thickness from equating the flow at 
the edge of the boundary layer to zero yields the equation 



When the value of 5 is substituted in equation (5.9) 


(5.10) 
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Computing the velocity of the flow of heat and returning to the 
previous variables yields in succession 


or 



(5.12) 


The heat given off by the sphere for small values of the Reynolds 
number S for air was considered by Leibenson. As a result of the 
accurate solution of equation (5.8), he obtained (reference 6, p. 280) 
equation (5.12) but with the numerical coefficient 0.85. 

S. free convection at plane vertical wall. - Let the origin of 
coordinates be at the edge of the plate; the x-axis is directed along 
the plane of the plate; the y-axis is perpendicular to the x-axis. 

From the physical considerations it is clear that the velocity and 
the temperature difference T - Tq (where T Q is the temperature of 

the air at large distances from the plate) differ appreciably from zero 
only in the thin boundary layer at the surface of the plate. 


For this reason, the process of free convection is described by 
the boundary -layer equations 


, Su’ 

u' v—p + v •, , 
ox' oy 1 


5u' _ c^u' 


Sy' 


+ g(3 (T' - Tq) 




u 


, ^ v , 


+ v ' Sr 


oj 


'Su* 3v' n 
X — r + 's — r = 0 

ox ' oy 1 

For the boundary conditions 


J 


( 6 . 1 ) 


u = 0 v = 0 

T = T i 

for 

o 

ii 

u = 0 

T = T 0 

for 

vs 

ii 

8 


( 6 . 2 ) 
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where is the temperature of the plate , 


From the parameters entering the equations, measures of length and 
■velocity can he constructed: 

r 11/3 1 


[eP < T i " T o>] 
fug^l - TqJpJ 1 ^ 


Equations (6.1) reduce to the form 


ox dy 


y OU , S 2 u o 

, S’ ^ - ry + * 

D oy 


•.* r Bu , i 

^jo s 


where 


d (C) = ] 

d (a.) = 0 


The boundary layer of convection 6 is next considered and the 
problem is solved by the approximate method. 


The first approximation is 


*- 1 "S 


7 

u = 2_ 
65 


y£ + 

2 3 


Substituting the values of and. b in the left side of the 
second equation of equations (6.4) yields the second approximation 
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The determination of the second approximation for the velocity, 
which is obtained without difficulty, is omitted. The boundary -layer 
thickness is determined from the condition (S^/oy ) g = 0. Thus 


R 3s 240 
5 5 = IIP 


or 


5 = 3.06 [ p 


1/4 


( 6 . 8 ) 


Substituting the value of 5 in equation (6.7) yields the temper- 
ature profile 


16 Ay 


* = 1 ‘ IT \5 ) T 11 \8 

whence the rate of heat flow will be 


i ri 


6 /y\ 5 10 /y> 4 


11 l 8 ) + 11 l 8 


(6.9) 


-(SL • “• (») 


1/4 


The accurate solution of the problem of the heat -transfer from a 
heated vertical wall was obtained by Pohlhausen (reference 7), who 
showed that if the magnitude 



( 6 . 10 ) 


is introduced as an independent variable and in place of the required 
functions, new functions defined by the relations 


u = 4UC 2 '\JT qp' (£) 


T - T 0 
*1 - T 0 


*U) 


are substituted, the initial system of equations takes the form 
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qp' • » +3 (pep" - 2<p’ 2 + £ 

+ 3P {p£ T = 0 

The solution of this system is obtained by numerical integration. 
From equations (6.10) and (6.11), the boundary-layer thickness will be 
of the order S ~ . In figure 4, taken from reference 8, is given 

the curve of d for air according to Pohlhausen. On the same figure 
the dashed curve shows the graph of d computed by equation (6.9); for 
y/5 <3/4 both curves practically coincide. 

7 . Stat ionary turbulent dif f ua ion . - In the solution of certain 
problems of geophysics, it is necessary to deal with the equation of 
diffusion in which the so-called- coefficient of turbulent exchange A 
enters instead of a constant coefficient of diffusion. For steady con- 
ditions, A depends only on the distance z. 


( 6 . 11 ) 


The problem is restricted to obtaining the second approximation, 
which is entirely sufficient for practical purposes. The problem of 
the character of the dependence of the coefficient of turbulent exchange 
on the distance z is not considered. Both sides of the diffusion 
equation 


u 




(7.1) 


are multiplied by A (z) and the new variable is introduced: 



Equation (7.1) then assumes the form 



(uA = 4> (q)) J 
Let the boundary and initial conditions be 

q» = q Q for z = 0 

q' = 0 for z = 00 

q’ = 0 for x = 0 


) 


(7.2) 


(7.3) 


(7.4) 
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In ord.er to simplify these conditions, q' = q^qj then 

q = 1 for z a 0 

At the edge of the boundary layer, the condition must exist that 
0, that is, that q (5) = 0. 


The first approximation has the form 

4 = 1 " 5 

The second approximation is given by 


= i - a + - 

5 5 


• r p 1 ! 

5 

5 3 

_Jo 


y(q - y)$ dy - ^ y(S - y) $ dy (7.6) 


21 = _ I + A 


S 5 2 


y $dy - - y(5 - y) <t> dy 


(7.7) 


The condition (dq/c5r])g = 0 gives for the boundary -layer thickness 
5 the equation 


$dy = S £ 


(7.8) 


From equation (7.7) for rj = 0 , 


^ .Air 

v°y n 5 s 3 

' ' 0 Jo 


y(6 - y)$ dy 


(7.9) 


By making use of equation (7.8), equations (7.9) and (7.6) for the con- 
centration may be transformed into the form 


■■■ hi ■■iiHimn ■ ■ iii mu !■ i n ■■ n i m lam n un i aw i 
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where 


r>8 





= I + _L 

^ O' 6 8 2 


0. 


d 1 

y$dy = ~L 
; 2 


r’l 


q. = 1 + 


y(ri - y) $>dy - tj 


(7.10) 


5 

J-i - I y 4> dy 


Jo 


J 2 = J y 2 $ dy 


Comparison with the exact solutions shows the effectiveness of the 
given approximation method. 


Translated by S. Reiss, 
National Advisory Committee 
for Aeronautics . 
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